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Abstract 

We use QCD sum rules to obtain the weak parity-violating pion-nucleon coupling 
constant f n NN- We find that f^NN ~ 3x 1CT 8 , up to an order of magnitude smaller than 
the "best estimates" based on quark models. This result follows from the cancellation 
between perturbative and nonperturbative QCD processes not found in quark models, 
but explicit in the QCD sum rule method. Our result is consistent with the experimental 
upper limit found from 18 F parity-violating measurements. 
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In this Letter, we use the method of QCD sum rules with the electroweak and QCD 
Lagrangians to predict the weak parity- violating (PV) pion-nucleon coupling constant, 
fnNN- The theoretical prediction of f n NN is an important and challenging problem. To- 
date, the most accurate PV experiments have only shown 1,2 ) that the upper limit for the 
magnitude of this coupling constant is 3-5 times smaller than the "best value" predicted 
by DDH 3 ) on the basis of a quark model and somewhat smaller than that in a similar 
calculation carried out more recently. 4 ) Since that time others have tried to estimate 
fwNN by means of chiral soliton models 5 ' 6 -* and QCD sum rules. 7 -* This coupling is of 
particular interest because of its sensitivity to the neutral currents contribution of weak 
nonleptonic processes at low energies. 2 ) 

QCD sum rules have been shown to be able to reproduce known properties of the 
nucleon, e.g., fx p , fj, n , gA, and of other hadrons. 8 ) However, they have rarely (if ever) been 
used to predict unknown properties. Keeping terms in the operator product expansion 
(OPE) up to dimension 5, we show that there are two main terms in the sum rule for 
fnNN- the unit operator and a dimension D=3 susceptibility. By using an analogous sum 
rule for the strong coupling constant, g^NN, to evaluate this susceptibility, we are able 
to determine the weak coupling f n NN- An important aspect of the present work is that 
we demonstrate that there is a cancellation between perturbative and nonperturbative 
QCD modifications of the weak process. 

We employ a two point function for the nucleon in an external pionic field. Our 
current is the usual one 9 ) 

Vp (x) = e^V^a^y^hV^*), 

fj P (y) = e abc [u b (y) lu Cu aT (y)}d c (y)Yl 5 , (1) 

where e abc is the antisymmetric tensor, C is the charge conjugation operator, and a, b, c 
are color indices. The neutron currents are similar, with the interchange of d <-> u. 
Since the most general weak PV n-N coupling is 2,3 ' 10 ) 

H PV (nNN) = ^p(r x 0^ , (2) 
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only charged pions can be emitted or absorbed. For definiteness, we consider the absorp- 
tion of a 7r + so that an initial neutron is converted to a proton, and the correlator we 
consider is 

n = ij d A xe tx - p < 0|T[77 P (x)77 n (0)]|0 >^+ . (3) 

The general form of II for the parity-violating pion-nucleon coupling, as dictated by 
relativistic invariance, is 

u PV = n e 1 + u p , (4) 

with p = 7^p M . 

The phenomenological evaluation of the correlator is carried out by deriving a dis- 
persion relation for II through the insertion of a complete set of physical intermediate 
states of spin | in the expression of Eq. 3. Using the usual terminology, we refer to this 
as the right-hand side (RHS). One finds for the parity- violating part of Eqs. 3,4: 

nPVf 2\rhs ^nUnn(p 2 + M 2 ) 

II; (p') ntls = — _ h continuum. (5) 

The double pole term corresponding to the insertion of the one-nucleon intermediate 
state in Eq. (3) has contributions both from the weak pion-nucleon vertex and the 
parity violation in the nucleon state itself. As will be shown below, in our microscopic 
calculation using the two-point form, only the latter process contributes. The parameter 
A at is related to the amplitude for finding three quarks in a nucleon at one point and has 
been determined in a number of sum-rule calculations. 8 - 1 M is the nucleon mass. As is 
usual in the method, the physical property of interest, f n NN, is obtained by treating the 
double-pole term explicitly, while the continuum and excited states are included in the 
numerical analysis via a parameterization, as discussed below. 

The microscopic evaluation of II, based on QCD and electroweak theory (the so- 
called left-hand sides (LHS) of the QCD sum rules for II), is obtained by means of a 
Wilson coefficient expansion in inverse powers of p 2 . In this work we keep diagrams up 
to dimension D = 5. The lowest dimensional diagrams which we consider are shown 
in Fig. 1. The higher dimensional diagrams which we include are obtained from those 



shown in Fig. 1 by the substitution of Figs. 2b and c for the pion-quark vertex, Fig. 2a. 
The propagators corresponding to the three diagrams of Fig. 2 are: 

oab %T ' 11 r , „,5 sab 

sf h = ^ T-Trg nqX n < m > $ ab i 5 , 

Sf c = -^rn-<qq>g nq r-nx 2 1 5 , (6) 

with XnQng < qq >=< 9*75? >tt and < q ij 5 g c a ■ Gq > w = tjIq < qq > r ■ n. Here g nq is 
the pion-quark coupling, which is not explicitly used in the present calculation, and G 
represents the gluon field. The susceptibility Xn enters in the evaluation of both strong 
and weak pion-nucleon coupling constants, while uIq enters only for the weak one. We 
will discuss the treatment of these parameters below. We only consider the even sum rule, 
namely that for Il e ; that for n o involves further unknown susceptibilities. The evaluation 
of the diagrams is straightforward. 

For the weak Hamiltonian, we take H w = ^| ( J M jJ + JVJVj) with 

N» = uY(A u + B ul 5 )u + dY(A d + B dl 5 )d, (7) 

where 9c is the Cabibbo angle and A u , Ad, B u , Bd, are given by 

1 8 
A u = -(1 - - sin 2 ^), 

A d = \ sin 2 M> 

B u = ~B d = — - , (8) 

with 6w the Weinberg angle. This is the standard model Hamiltonian, which we use for 
the main part of the calculation. We then discuss the QCD effects on our results. 

Since momentum can be transferred in the weak point-like interaction, shown by 
wavey lines representing Z° in the figures, there is an additional integral to be carried 
out in the evaluation of n. For example, we obtain for Fig. la 

n- = i| a**. J { ^m±^ (AA _ AMy^d'y , (9) 
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There is no PV contribution from Figs, (lc) and (Id), and the sum of Figs, (lb) and 
(le) vanish. The integrals in Eq. (9) are evaluated by standard Feynman techniques, 
with dimensional regularization. The result is 



fla . 2 -12 G F g„ qK _ fe^d^x ( 2 , n _ 22 



n eV) = limy {- - + 2 7 ~ y ) IV?„ - -4,,^) 

Gf /22 2\ g . 2 



^ — - 2 7 + - )p b ln( V) - A u 5 d ) , (10) 



3 • 2% 6 ^ 9 V 3 

where the infinite term oc e _1 is assumed to be regularized away, and 7 is the Euler 
constant. Taking the Borel transform one obtains 

= ^^(»/3-7) «S. (ID 

where M B is the Borel mass. The other diagrams can be evaluated in the same manner. 
The results from the processes of Figs. 1, and those obtained from Figs. 1 with the 
substitution of Figs. 2b and c for Fig. 2a are 

12tt 2 



+ \x* aL~^ 9 M 2 E 2 -ml a E X L~^\ 



M 2 



= U NN \ 2 N e- M2 l M H2 WB -l) (12) 

where a = — (2n) 2 < qq > and = (2ir) A }? N / g nq . The lowest dimension pion-quark 
vertex renormalization diagrams for f n NN are shown in Fig. 3a. In our approximation of 
a contact weak interaction, the contribution of Figs. 3a-c vanish under a Borel transfor- 
mation. The mechanism of Figs. 3d and e do not appear in the external field method. We 
do not include gluon condensate diagrams for J-^nn] they are of the same order or smaller 
than uncertainties of our calculation. The factors containing L, L = 0.621 ln(10M B ), give 
the evolution in Q 2 arising from the anamolous dimensions, and the Ei(M^) functions 
take into account excited states to ensure the proper large-M^ behavior. The last line 
in Eq. (12) is the Borel transform of the double-pole term from the phenomenological 
(right-hand) side, Eq. 5. The direct proportionality to sin 2 9 W should be noted. 
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Finally, by explicit calculation or Fierz reordering, we can show that the contribu- 
tion for W ± exchanges vanish. Thus, as required by symmetries 3,10 ^, we find no charged 
current contribution to the weak PV pion-nucleon vertex; such a contribution requires 
strangeness-changing currents and would thus be reduced by sin 2 6 C ~ 0.05. Since we 
neglect strangeness in the nucleon and strangeness-changing currents, we obtain no con- 
tribution. 

As we shall demonstrate below, the first two terms in the theoretical form for IT e 
given in Eq. 12 are of opposite sign and tend to cancel. This is a crucial point. For this 
reason it is essential to either determine the value of the susceptibility x-n from g w NN or 
to eliminate it from our equations. We do both as an aid in determining the stability 
of our solutions. First, we determine x-n directly in terms of g^NN [as a function of the 
Borel mass] by using the sum rule for the strong coupling, which is analogous to Eq. 12, 
and attempt to use the result to determine f n NN- Second, we eliminate Xtt from the PV 
and strong coupling sum rules and find that we can determine f^NN in terms of g n NN- 
Details are given below. 

We use the correlator given by the two-point function of Eq. 3 for the strong as well 
as the weak interaction. The general form differs from Eq. 4 by the presence of a 7 5 
in each term. The phenomenological (RHS) for the strong pion-nucleon coupling is now 
given by 

K{p 2 )RH s = ^f_^ 2 7 5 + continuum. (13) 

Unlike the weak PV pion-nucleon coupling, the evaluation of the strong one leads to 
a problem in that there is no double pole on the right-hand (dispersion relation) side. 
However, as shown by Reinders et. al. 11 ), the value of the coupling constant g^NN found 
in this way is virtually the same as that found by means of a 3-point function, which 
circumvents the lack of a double pole problem. 

Keeping terms up to D=6, shown in Fig. 4, for the theoretical side (LHS), and taking 
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the Borel transform we obtain the sum rule for the strong pion-nucleon coupling: 

9«nn\ 2 n e" M2/M - = M%L'^E 2 - M^aL^E, + \a 2 L^ 

+ < 9l ° 2 g > E * Ml - < ft? > E MUf - lnM% (14) 

where < g 2 G 2 > is the gluonic condensate. 

Before we discuss our detailed evaluation of the sum rules to obtain our estimate of 
fwNNi let us discuss the structure of Eqs. (12) and (14). First, as we discuss below, if 
we use PCAC to evaluate x-k we find that \- K (i = - 44 GeV 2 . With this value, the x-k term 
dominates both Eqs. (12) and (14), with the result that g n NN — 155 [in contrast to the 
experimental value of 13.5]. With this value of x-k we find that f^NN > 1CT 6 , at least 
an order of magnitude larger than experiment. 

Secondly, since Xn is the only unknown in Eq. (14), we can estimate the vacuum 
susceptibility using the experimental value of g n NN = 13.5: this gives x-k a — -0.94 GeV 2 , 
two orders of magnitude smaller that the PCAC value [see discussion below]. With this 
value one finds that the first two terms in Eq. (12), the leading terms for f n NN, almost 
cancel. Note that the second term involving x-k enters with the opposite sign in the two 
equations for f w NN and g n NN, respectively. This is the source of the very small parity- 
violating pion-nucleon coupling in comparison with quark model: there is a cancellation 
between the dimension zero model-like term using perturbative quark propagators and 
the vacuum pion susceptibility term. 

However, we find that the sum rule obtained for f^NN [Eq. (12)], using the value of 
X 7r (M 2 3 ) extracted from Eq. (14), is not stable in M^. Therefore we cannot obtain a 
reliable estimate of f n NN by this method. 

We find that we can obtain a satisfactory sum rule to determine f n NN by eliminating 

X-k from both Eq. 12 and Eq. 14 by taking derivatives with respect to M|. With this 

procedure, and taking the ratio of the weak to the strong sum rule we obtain the new 

sum rule for the weak in terms of the strong coupling constant: 

Unn 2 [Ml - AM%){EkM% - amlE,) 

g nNN w N (2M* N + 3M4 -9M 2 N Mi)(l2E 2 M 4 B + 3<G*> E o y { ) 
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where c w = G F sin 2 6 w (f - 7)/(12tt 2 ) = 6.7 x 10~ 8 GeV~ 2 . The sum rule is quite 
stable with a plateau in M 2 B in the region expected, as shown in Fig. 5. Because of the 
strong cancellation between the first two terms in Eq. (12) [dimension and dimension 
2 terms], the dimension four term with the unknown parameter nig is important for the 
final numerical value of f w NN- We have taken mj = in Fig 5. Guided by the value of 
the parameter m needed in the nucleon sum rule 8 - 1 , we evaluate the sum rule given in 
Eq. 15 with mj taken over the range 0.0 to +0.8 GeV. From this procedure we find: 

Unn « (2.3 to 5.8) x 10' 8 for 
ml = (0to0.8)GeV. (16) 

For negative values of mj the value of f n NN becomes smaller and even negative, but we 
did not find stable solutions for sizable negative values of this unknown parameter. To 
be consistent with the neglect of gluon condensate terms we quote as our central value 
of f-wNN that with = 0, shown in Fig. 5, as 

f wNN » 2.3 x 10- 8 . (17) 

This coupling constant is an order of magnitude smaller than the "best values" of Refs. 3 
and 4. As emphasized earlier, this result follows from the cancellation of the two leading 
terms in Eq. 12. The first LHS term in that equation, a unit dimension term which 
would correspond to a quark model type calculation gives, a value for f^NN ~ 2 x 10~ 7 , 
similar to the quark model value. The second term, involving the nonperturbative QCD 
vacuum susceptibility, Xn, strongly cancels the first term. Because of this cancellation, 
we cannot expect Eq. 17 to be very accurate, but we find a clear explanation for the 
small value of f-nNN, consistent with experiment. 1 ' 2 -* 

The results given in Eqs. 16 and 17 have been obtained using the Hamiltonian of 
the standard model (see Eqs. 7 and 8). Let us now consider the strong interaction 
modifications. These have been estimated in Refs. 3 and 4 using the renormalization 
group method. In the notation of Ref. 3, the operators involved in our calculation are O4 
and O5. Using the tables in Refs. 3 and 4 we find that the our parameter for the parity 
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violation, A&B U — A u Bd, would be changed by less than a factor of two in magnitude. 
Since the same parameter appears in all terms, this gives the overall uncertainty arising 
from strong interaction modifications. Therefore, the main conclusion of our work is not 
changed. 

There are two relevant features that we would like to point out. The first one is that 
the use of pseudovector coupling also circumvents the problem of a lack of double pole 
for the strong interaction constant. For the Lagrangian 

C„ NN = ^™ $ N tYl 5 r • ^V M 4 (18) 
m n 

we can treat V M 0,r as a constant external axial vector field. The QCD sum rule is then 
identical to our calculation of qa- 8 ^ At the quark level, we have 

L nqq = ^ T i' q h^ 5 rij q Vj n , (19) 
where /„- is the pion decay constant. From our previous result for ^ we then obtain 



9ttNN 9a 

, 2M g A M 

9nNN = 9-kNN = —f (20) 

which is just the Goldberger-Treiman relation. 

As a second feature we wish to attempt an independent estimate of Xn- For this 
purpose we first use PCAC to obtain 

_ — f rri^ 

< 0\u i^u - d i lb d\it° >= " - e~ iq - x . (21) 

We then use the work of Belyaev and Kogan 12 ), which assumes saturation of a sum by 
one pion states: 

< 0|gi7 5 r 3 g|0 > 7r = -^L <j> n J d 4 xe lQx < 0\ui>y 5 u - di^d\n >< ir\q i>y 5 T 3 q\0 > Q ^ 
= ^ %f = xA.<W> (22) 



9 



As described above, the value of Xn obtained in this manner is more than an order 
of magnitude larger than that found by using the value of g n NN from experiment. Once 
more we point out that if we use it in Eq. 14 we find an order of magnitude discrepancy 
with the strong coupling constant, g n NN- Furthermore, it is clear that this value of 
X-k is inconsistent with Eq. (12), since by eliminating it with derivatives with respect 
to the Borel mass we obtain results an order of magnitude different than with its use. 
We conclude that Eq. (21) cannot be correct. We are not certain where the method of 
Belyaev and Kogan errs, but we believe that it is suspect. 

In conclusion, we find that the weak PV pion-nucleon coupling due to neutral cur- 
rents is as small as that due to charged currents, ~ 3 x 1CT 8 . This result agrees with the 
conclusion of the chiral soliton model of Kaiser and Meissner 5 ^ , but not that of Kaplan 
and Savage 6 ^. Our result also disagrees with quark model calculations 3 ' 4 ) and with a 
previous QCD sum rule calculation. 7 ) If the coupling is as small as we estimate, it cannot 
be separated from the charged current contribution and thus cannot be found experimen- 
tally. Although we have omitted gluon condensate corrections to the PV correlator, our 
result is sufficiently small that these corrections will not alter our conclusion. Finally, we 
point out that in the two-point QCD sum rule method used here, the small value of f n NN 
which we obtained is the result of a cancellation between a process which can be treated 
in quark models and a vacuum process identified in the method of QCD sum rules. 
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Figure Captions 

Fig. 1. Lowest dimension quark diagrams for the PV weak pion-nucleon vertex. The 
dashed line represents a charged pion and the wavy line a Z°. 

Fig. 2. Quark propagator modifications in an external pion field. 

Fig. 3. Pion-nucleon weak vertex correction diagrams. 

Fig. 4. Diagrams contributing to the calculation of g n NN- 

Fig. 5. Solution for f n NN from Eq. 15 as a function of M B . 
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